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ABSTRACT 

I decompose the unstable growing modes of stellar disks to their Fourier components and present the 
physical mechanism of instabilities in the context of resonances. When the equilibrium distribution 
function is a non-uniform function of the orbital angular momentum, the capture of stars into the 
corotation resonance imbalances the disk angular momentum and triggers growing bar and spiral 
modes. The stellar disk can then recover its angular momentum balance through the response of 
non-resonant stars. I carry out a complete analysis of orbital structure corresponding to each Fourier 
component in the radial angle, and present a mathematical condition for the occurrence of van Kampen 
modes, which constitute a continuous family. I discuss on the discreteness and allowable pattern speeds 
of unstable modes and argue that the mode growth is saturated due to the resonance overlapping 
mechanism. An individually growing mode can also be suppressed if the corotation and inner Lindblad 
resonances coexist and compete to capture a group of stars. Based on this mechanism, I show that 
self-consistent scale-free disks with a sufHcient distribution of non-circular orbits should be stable 
under perturbations of angular wavenumber m > 1 . I also derive a criterion for the stability of stellar 
disks against non-axisymmetric excitations. 

Subject headings: stellar dynamics, instabilities, methods: analytical, galaxies: kinematics and dy- 
namics, galaxies: spiral, galaxies: structure 



1. INTRODUCTION 



Both iV-body simulations (jHohl I Il971f ) and analyti- 
cal methods (Kalnajs 1978; Jalali & Hunter 2005, here- 
after JH) show that global instabilities can generate 
barred structures in galactic disks. Apart from the bar 
mode, which is an isolated event in frequency space, 
global spiral modes seen in the eigenspectra of cored 
stellar disks (Jalali 2007, hereafter Paper I) consti- 
tute a discrete fam ily that bifurcates form stationary 
Ivan Kampen I (|1955l ) modes. But not all spiral struc- 
tures are g lobal modes as disturba nces induced by close 
neighbors (Bvrd fc Howard II 1992f ) and density inhomo- 
geneities tToomre II1990D inav also create the spiral pat- 
terns of the observed galaxies. 

A mode of a stellar disk is a mathematical entity that 
comes out of an eigenvalue problem. However, its phys- 
ical origin in isolated systems has not yet been under- 
stood clearly. Lynden-Bell & Kalnajs (1972, hereafter 
LBK) attempted to explain a mode through the trans- 
port of angular momentum between different parts of 
the disk. They suggested that the inner Lindblad reso- 
nance (ILR) releases the angular momentum of central 
regions and the spiral structure transports it to the outer 
parts through the corotation (CR) and outer Lindblad 
resonances (OLR). This mechanism is favored by some 
galactic dynamicists (,Athanassoula ,,2003, ), but it is se- 
riously challenged by Toomre's (1981) theory that says 
that feedback through the galactic center is a critical in- 
gredient for growing modes. In Toomre's theory, on the 
other hand, the modeling of feedback as the reflection of 
a leading spiral wave at the galactic center and its emer- 
gence as a trailing one, is a simple description of a very 
complicated dynamics that governs the motions of stars. 
Unresolved issues concerning the evolution of unstable 
modes include the following: (i) The swing amplification 



theory is not capable of predicting the fate of a growing 
mode against other stationary and unstable modes that 
coexist in the eigenspectrum of a given model, (ii) How 
does the bar mode saturate? (|e.g.. Khoperskov et al. I 
|2007| ). (iii) Why does the nonlinear bar terminate al- 
most at the corotation radius? (Sellwood 1981) (iv) We 
should also understand the origin of different species in 
an eigenspectrum and interpret their continuous or dis- 
crete nature, and the distribution of their pattern speeds 
and growth rates. 

Stellar orbits in a galactic disk begin to evolve once 
the surface density deviates from its equilibrium state, 
'Eq(x), and develops a time-varying mean-field potential 
Vi{x,t). Here x denotes the position vector of stars at 
the time t. In the linear regime, we are usually interested 
in density waves that grow/decay according to the expo- 
nential law e''* and rotate with the fixed pattern speed 
Op. For a wave of m-fold symmetry, the perturbed po- 
tential becomes 



Vi = ee'*V {x^mflpt) . 



(1) 



Since we are dealing with infinitesimal perturbations, I 
have introduced the small parameter e so that ee^* ^ 1. 
From II]) one arrives at the equations of motion 



x = V, 



V = — 



dx 



dx ' 



(2) 



where Vo{x) is the equilibrium potential field generated 
by galactic stars and a possible dark matter halo. In 
writing equation ([2]), I have assumed that the motion 
of stars is restricted to the disk plane. When the equi- 
librium state is axisymmetric and the dark component 
is spherical, Vq becomes a function of radial distance to 
the galactic center and the unperturbed equations (with 
e = 0) are integrable. In such a circumstance, the phase 
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space is filled by rosette orbits denoted by [xo{t), Vo{t)]. 
The growth of perturbations, whatever the magnitude of 
ee*** <C 1 may be, deforms stellar orbits. Orbital defor- 
mations are measured hy x = x — Xq and v = v — Vq, 
which can be used in (l2l) to obtain 



dx 



dv 





■ i - ee"' 


dV 


dx 


dx^ 





(3) 



Although a proper equilibrium distribution function 
fo{x,v) can self-consistently reproduce So(a;) using 
rosette orbits, the perturbed density T,i{x,t) (corre- 
sponding to Vi) cannot be supported by rosette orbits 
alone and orbital deformations are necessary for the self- 
consistency of density waves. According to equations ([2]) 
and ([3]), orbital deformations of O (ee"*) are sufficient to 
support the growth of density/potential perturbations 
up to the same order of magnitude of such deformations 
over a time scale of I/O (ee"*). As the time is elapsed, 
the amplitude of perturbations increases exponentially 
and the solution of the linearized coUisionless Boltzmann 
equation (CBE) fails when ee** ^ 1. During my mode 
calculations, I realized that the orbital axes of certain 
stars librate in a coordinate frame that rotates with the 
density pattern. This resonant capture initially seemed 
to be a higher-order nonlinear effect but further exper- 
iments showed that the resonant gap is constrained by 
the magnitude of density perturbations. The complex 
behavior of stars for infinitesimally small yet non-zero 
ee** <C 1, and the role of resonant stars in the generation 
of discrete galactic modes, are investigated in this paper. 

I use the results of Paper I and introduce a new dynam- 
ical mechanism that sparks unstable mod e s and governs 
the singular oscillations of Ivan Kampen I (|1955l ) modes. 
Resolving the origin of instabilities and amplitude sat- 
uration precede my nonlinear calculations, which were 
made feasible in Paper I by the Petrov-Galerkin method 
and reducing the CBE to a system of nonlinear ordinary 
differential equations. Those reduced equations, how- 
ever, are valid only when orbits are regular and averaging 
over angle variables is allowed. As unstable modes grow, 
chaotic orbits come into existence and the weighted resid- 
ual form of the CBE must be modified to handle them. 
I quote some of the results of such modifications in this 
paper when I discuss the issue of mode saturation. In 
Paper III, I will give a full account of the mathemati- 
cal and numerical modeling of stochastic layers, and will 
analyze modal interactions after their saturation phase. 

For the cored exponential disk embedded in the field 
of the cored logarithmic potential, I describe the de- 
composed Fourier components of unstable bar and spiral 
modes in S|3] and highlight the existence of a phase shift 
between different components. In 21 1 derive a condition 
for the corotation of the orbital axes of an ensemble of 
stars and explain the role of such a synchronous motion 
in pattern formation. I dedicate Sj5] to exploring the or- 
bital structure of a perturbed stellar disk and identify a 
resonance mechanism that can generate both stationary 
and growing modes. I reveal the mechanism of angu- 
lar momentum transfer between Fourier components and 
derive analytical expressions for the growth of resonance 
zone. I address the origin of instabilities in fj6l present 
a saturation mechanism for unstable modes in ^ and 
discuss about the global stability of soft-centered and 



scale-free disks. I explain the restrictions of LBK's mode 
mechanism in ^J8] and end up the paper with concluding 
remarks. 

2. THE MODEL 

The calculations of the present study are carried out for 
the cored exponential disks of JH whose eigenfrequency 
spectra and mode shapes have been completely explored 
in Paper I. The model has a dark matter halo and the 
motion of stars in the equilibrium state is governed by 
the cored logarithmic potential 



Vo{R) 



:lnA/l 



(4) 



which is the resultant gravitational potential of luminous 
and dark components. Throughout the paper, all length 
and velocity variables are normalized, respectively, to 
the core radius of the potential {Rc) and the asymp- 
totic velocity of stars on circular orbits (vq) by setting 
Rc = ''^o = 1- To describe the physical quantities in the 
configuration space, I will use the usual polar coordinates 
{R, (/)) and their Cartesian counterparts {x, y) where R is 
the radial distance from the galactic center and (j) is tbe 
azimuthal angle. The specific model that I adopt here is 
a relatively cold, near-maximal disk with no dark matter 
concentration in the region with rising rotation curve. 
The model parameters are set to (iV, A, a) = (6, 1,0.42) 
where N is an integer exponent that controls the propor- 
tion of circular orbits and the disk temperature. Larger 
values of N give rise to colder disks. For = 6, the 
parameter Q of Toomre (1964) is marginally larger than 
1. The parameter A is defined as the ratio Rc / Rd with 
R]j being the core radius of the equilibrium density, a is 
a factor that controls the total mass of the stellar compo- 
nent. The chosen model with a = 0.42 is near-maximal 
within a radius of « 2.5Rc and dominated by dark mat- 
ter beyond it. The calculations of Paper I revealed the 
eigenspectrum of this model that includes a compact bar 
mode Bl and a sequence of spiral modes SI, • • •, S6. I 
select modes Bl and S2 as my case studies of sections [3] 
through m The reason for choosing mode S2 is its exten- 
sive and prominent spiral structure and modest growth 
rate. In section [71 I will also display some results for 
modes SI, S3 and S6. 

3. DECOMPOSITION OF UNSTABLE MODES 

One of the advantages of the method developed in Pa- 
per I is that the density function of a mode can be read- 
ily decomposed to its constituent Fourier components in 
angle-action space. The phase shifts between density 
components determine the magnitude and direction of 
the torque that is exerted on each component. So we 
can probe the transfer of angular momentum and iden- 
tify the direction of its flow once a global mode develops. 

The perturbed distribution and Hamiltonian functions 
of an unstable mode can be expanded as (Paper I) 



C30 

/i(e,j,<)=Rc^ 

m,/ — — oo j 



oo oo 



Hi{e,J, t)=Re^ ^^bf ( J)e'('"^*+'^«) , (6) 

m,^ — — oo J— 

where 9 = [9r, 9^) and J = {Jr,J^) are the angle and 
action variables, respectively. The angles are defined 
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based on the radial and azimuthal frequencies 
of stars on rosette orbits so that 

0R = nR{j), d^^n^j). (8) 

A dot stands for the time derivative and TLq{J) is the 
integrable Hamiltonian of the axisymmetric equiUbrium 
state. For a normal mode the amplitude functions 

and depend on the time variable t through the 

simple exponential law exp(— iwt) with w = m^lp + is 
and I = V— 1. Here Vtp and s are the pattern speed and 
growth rate of an unstable mode of angular wavenumber 
m. On the other hand, the perturbed potential function 
Vi and its associated surface density Si can be expanded 
in the configuration space as 



Vi(i?,</>,<)=Rc^ ^ea™(t)Vrl(i?)e 



(9) 



a— — oo 

where crj™'(i?) and are surface density and po- 

tential basis functions, respectively. They satisfy Pois- 
son's integral and the bi-orthogonality condition 

oo 

D,{m)5j,j,5ra.,n,' = 2^ j 4,^p\R)a^p'\R)RdR. (11) 



Here Sm,m' is the Kronecker delta and Dj{m) are some 
constants that depend on our special choice of basis func- 
tions. Noting Vi =Hi, one can equate ^ and dH), mul- 
tiply the identity by exp[~2(Z0fl + m9(j,)] and integrate 
over the angle variables to obtain 

Y^bfmfiJ) = ^af (J), (12) 



TV 

*f (J) = i J ^P^"'^ (R) cos [WR + m {6^ -<f>)]deR. (13) 



Setting *J''(J) *™'(J) gives bf{t) = af{t), which 
is a remarkable simplification. A relation between a^{t) 
and d™' (t) then follows from the weighted residual form 
of the fundamental equation 

fi{Q,J,t)dJdQ = '>:i{R,(j),t)RdRd(j). (14) 



I obtain 



2 °° o° 

<('^^:d>^ ^ E^trff (15) 



l=~oc k=0 



(16) 



where 



Let C be the total angular momentum of the disk, 
which must be conserved in an isolated galaxy. This 
means that the torque dC/dt must vanish. Defining 



df{t) = e 
one obtains 

dC 



3 ' J V ' 
lujt (jjrnl I ^-\/-'nil\ 



dt 



^2rrm^e^e^'' ^ L„(0 = 0, 



j,fc=0 



(18) 
(19) 

(20) 
(21) 



The constant vectors (C/r'>^™0 and {uj',v^) are ob- 
tained from the linear eigenvalue equations of Paper I 
and a subsequent use of equation (fT5|) . Equation (|20p is 
analogous to equation (B6) in JH and Lm{l) shows the 
rate of angular momentum flow into/from the ^th Fourier 
component (in the radial angle 9r) whose corresponding 
density in the configuration space is 



oo oo 

Ei(i?,</>,t)=Rc^ ^ea7(^)al.'"l(i^)e™^ (10) E'i 



47r^ 



:ee 



j,k=0 "'^ ^ 



A 



V"'(i?)d^'e''"('^-"''*^ (22) 



The angular momentum conservation of the disk implies 
that some Lm{l) take negative and some other positive 
values with a total vanishing sum. This means that Yj\^ 
the pattern corresponding to Lm(l), is subject to a pos- 
itive torque from other components if Lm{l) > 0, and a 
negative torque otherwise. 

Figure [1] shows the variation of ^2(0 versus I for modes 
Bl and S2. The results are in accordance with the bar 
charts of JH: components with I < lose angular mo- 
mentum and those with I > gain it. As JH had already 
pointed out, a few components ensure the convergence of 
Fourier series in the 0ij-direction. The superposition of 
the I > components of the perturbed density, defined 



as 



4+(i?,0,i) =X]sUi?,0,t), 



(23) 



will thus be a pattern that experiences a positive torque 
exerted by 

-1 

J:[-{RA,t) = (24) 

/ — — oo 

In response, the reaction torque of "E^^ will drain the an- 
gular momentum of S^" so that C is conserved. Figure [2] 
shows modes Bl and S2 and their components that emit 
and absorb angular momentum. It is seen that in both 
modes, the phase of lags that of by a magnitude 
of « 90° and the positive parts of S']^~ fill the regions 
of negative E^^. Due to this phase shift, the angular 
momentum is transferred between I < and I > com- 
ponents and a counterclockwise torque is exerted on E'^^ . 
This phenomenon is more obvious in mode Bl that has 
a definite spatial orientation. The most important ques- 
tion relevant to the origin of instabilities arises now: Are 
the density components Y^[{R,(t),t) generated by stellar 
orbits at resonances? Analytical calculations of sections 
|4] and [5] provide the answer. 



4 



M. A. Jalali 




Fig. 1. — The angular momentum content of different Fourier components of modes Bl and S2. The sum y^."t_°° 



Lm{l) vanishes because 



the total angular momentum of the disk is conserved. The function Lm (l) is normalized so that the sum of positive components is unity. 
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Fig. 2. — Modes Bl and S2 and their components that release and absorb angular momentum. Left panels show the mode shapes, which 
rotate counter-clockwise. The pattern speeds and growth rates of modes Bl and S2 are (Op, s) = (0.918, 1.160) and (Hp, s) = (0.454, 0.216), 

respectively. Middle and right panels demonstrate the Sj~ and Ej"*" componenets, respectively. Solid lines in top panels highlight the 
orientations of the bar mode and its components. The prominent phase shift between E'^~ and E'j^^ is responsible for the gravitational 
torque between them. The isocontours show the positive part of the density from 10% to 90% of the maximum with increments of 10%. 



4. SYNCHRONOUS PRECESSION OF THE ORBITAL AXES 

The perturbed density Ei of a mode with the pattern 
speed Op must be supported by the slow motion of stars 
in a rotating coordinate frame of angular velocity flp. In 
other words, by sitting on a moving coordinate system 
one must be able to identify a group of stars that are 
losing their angular velocity, and another group of the 
same kind but with different initial conditions, gaining 
it. A density peak is expected in the region that the an- 
gular velocity of the slow ensemble is minimum. Since 
the angular velocity of a star is minimum at its orbital 
apocenter, the precession rate of the orbital axes of par- 



ticipating stars in the pattern formation should be in 
resonance with the pattern speed. The above scenario is 
legitimate as long as pattern stars stay far from Lindblad 
resonances that can generate higher density regions due 
to the radial slowing of stars. My computations show 
that all unstable modes fulfill this requirement and az- 
imuthal slowing is the main origin of density perturba- 
tions (see SjS]). 

The evolution of (j){t) is crucial for understanding the 
slow dynamics in the configuration space, but we usually 
have the perturbed distribution and potential functions 
in the angle-action space. It is thus useful to find how ((> 
depends on (8, J). The simplest relation can be obtained 
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by expanding exp[zm(0 — 9^)] in Fourier series of Or as 



+ OC 



where 



(25) 



(26) 



The functions S,mi{J) are computed by integrating along 
rosette orbits (over a period of radial oscillation), and 
they vanish for circular orbits when / ^ and for ra- 
dial orbits when m ^ —21. Multiplying equation (I25p 
by exp[— zm((/) — 0^)] and integrating the identity over a 
period of 9ji yields the following useful relation 



-|-oo 



(27) 



^— — oo 



Moreover, taking the partial derivative of (|25|) with re- 
spect to 9fi, multiplying both sides of the resulting equa- 
tion by exTp[—tm{(f) — 9^)] followed by an integration over 
9r, leads to (Scott Tremaine, private communication) 



-t-oo 



(28) 



^— — oo 



I now multiply (I25p by exp(im00), differentiate both 
sides of the resulting equation with respect to t and ob- 
tain 



-ReE 



+ 00 



dJR 



dJd, 



Re 



E Cm/ ( 



19, 



m9, 



(29) 



The perturbed motions of stars are governed by Hamil 
ton's equations 

dUi ■ dUi 



e = n(j) 



dJ 



that can be used to rewrite 



in the form 



(30) 



4> = Re E [l^R + ™^ 



-ReE 



/ — — OO 



4,\e 



(31) 



where [•••,•••] denotes a Poisson bracket taken over the 
angle-action space. On substituting for the complex con- 
jugate of exp(im(/)) from (|^5|) in (|3ip and noting that the 
second term on the right hand side of (pij) is of O (ee***), 
one finds 



l.k— — oo 

(32) 

which shows the functional dependence of </> on the ac- 
tions and the radial angle. 
Equation ([5^ can be written as 

m4,^mVl^{J)+f]{J,9R)+0{ee'') , (33) 



where the constant part mfi^ is obtained by setting I — k 
in the double summation of ([5^ and subsequent appli- 
cation of relations (P7)) and ([55)1 . r]{J,9R) is a periodic 
function of 9r that stands for all other terms of (|32p 
with I ^ fc. ri0(J) is thus the precession rate of the 
orbital axis of stars whose energy and angular momen- 
tum correspond to the action vector J. Those stars can 
contribute to a developing density peak when they linger 
at the apocenter of their orbits. Such a mass deposition 
will continue in a rotating frame of angular velocity Qp 
and over the time scale 1/0(66**) if the condition for 
synchronous precession 



holds. By defining 



(34) 



(35) 



^i„liJ) = l^niJ) + mVL^{J) - rnQp, 

and using p3|) . the condition (|34p takes the convenient 
form 

M°,(J)«0(ee**), (36) 

which is the definition of the CR. The ILR and OLR oc- 
cur if l/ijj^l diminishes for I = —1 and ^ = 1, respectively. 

A graphical representation of ([36]) will help to sharpen 
our understanding of the ensemble of stars that can sup- 
port a rotating pattern. Since there is a one-to-one, onto 
and invertible map J ^ for initially axisymmctric 
disks (excluding Keplerian and harmonic oscillator po- 
tentials), any functional form of the actions can be de- 
scribed in terms of $7 as well. For instance, one may write 
/z^(J) = /i55^(f2). I therefore identify resonant regions in 
the frequency space. For demonstrating the distribution 
of dependent quantities, I will use the pair {^r, f2i) as the 
coordinates where fli — fl^ — ■^^r- The reason for this 
choice is that the orbital frequencies of soft-centered stel- 
lar disks fill a very narrow region in the fl-space (Hunter 
2002, hereafter H02; Figure 1 in JH), which does not 
provide enough resolution for the visual identification 
of some fine structures. Denoting flo — [f^^Jmax ^ 
cored stellar disk, the center of the disk and infinity cor- 
respond to (nR^n,) = (2f]o,0) and inR,ni) = (0,0), 
respectively. Figure [3] shows the frequency space of the 
cored logarithmic potential. The lower boundary (f2fl- 
axis) corresponds to radial orbits with = and the 
upper curved boundary fli — Tc{i^R) is determined by 
the orbital frequencies of circular orbits with Jr = 0. 

I have plotted in Figure[3]the countours of |/i°j(J7)| for 
m — 2 and Hp — 0.454ilo, which is the pattern speed of 
mode S2. Lighter regions mark the stars that are closer 
to the exact CR. The intersection of the straight line 
^2(^^) = with the frequency space corresponds to the 
stars at the CR. The centerline of the CR zone begins 
at the location of a star moving on the corotation circle, 
and extends into inner regions. 

By decreasing Hp, the CR zone is pushed to the out- 
skirts of the disk where orbital frequencies are small and 
the surface density has dropped substantially. Further- 
more, the line /z^(r2) = will not intersect the frequency 
space, and there will not be a CR if Hp > Hq. For a sim- 
ilar reason, the ILR will be absent if 

(37) 

Hr . (38) 



r / in \ 
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Fig. 3. — The frequency space of the cored logarithmic poten- 
tial with Rc = VQ = Qo = 1- The upper boundary Qi = rc{f^R) 
corresponds to circular orbits with Jjj = 0. Contours show the 
variation of the CR indicator |^2(n)| over the frequency space 
for Qp = 0.454. Lighter regions correspond to smaller values of 
l/ijC^)!- The function ^§(f2) = 2{Q^ — Qp) vanishes at the inter- 
section of the drawn straight line and the frequency space. The star 
at the intersection of the line /ijCfJ) = and the upper boundary 
moves on the ixirntation circle. _ 

Inequality (|37j) is satisfied when the hne /i~ (H) = 
does not cross the boundary curve fli — Tc{i^B)- These 
properties appear to correlate with the computed eigen- 
frequencies of unstable modes: according to the eigen- 
spectra of Paper I, the pattern speeds of growing modes 
lie in the interval 



(39) 



The lower limit may be violated by cold disks and disks 
with inner cutouts (see i j7.2p . The upper limit is vio- 
lated by some rapidly rotating modes in galaxies with 
dark matter, and by mode Bl of cutout disks. After ex- 
ploring the orbital structure of perturbed disks in SjSj I 
will explain the constraint (j39p and the exceptional cases 
that may not satisfy it. 

5. STELLAR DYNAMICS IN THE PERTURBED DISK 

The orbital dynamics of stars can be understood 
through investigating the induced dynamics by each 
Fourier component. We usually split the phase space 
to several subspaces by using a Fourier expansion in 9 a, 
and each subspace is associated with a pattern compo- 
nent (e.g., Figure [2|). In the linear regime, the perturbed 
motion of stars due to the Ith component can be traced 
by carrying out a canonical transformation 

i0R,O^,jR,J4,) ^ {wi,W2,Il,l2) , (40) 

defined by the generating function ()Lvnden-Bell |[l993 ) 
S ^ {WR + m0^-mnpt)h+eRl2. (41) 
This yields the transformation rules 
9S , dS 

and the new Hamiltonian 

IC = Ho{I)+Hiiw,I,t) + dS/dt, (43) 

so that 



ic=no{i)-mapii 



^^mk^j-^^i[wi+{l-k)w2] ^ (44) 



k— — oo j—Q 

From the equations of motion 
dJC . 
ow 



dK. 

a/' 



(45) 



one can verify that W2 — ^rI + 0{ee^*) is an increasing 
function of time, which can be made a cyclic coordinate 
by averaging /C over W2- Consequently, the action I2 = 
Jr — III = Jr — IJip/rn becomes an adiabatic invariant 
so that I2 = 0(6^6^"*), and the dynamics is reduced to 
the flows governed by the averaged Hamiltonian 

^(m, l)={iq^,^lCo + ee^Ti , (46) 

lCo{I)=no{I)~mnph (47) 

lCi=AUl)cos[wi+dUl)], (48) 

where I have kept the real (physical) part of {)C)w2 ^^J^d 

<,(/) - VX^ + Y\ dUl) = arctan , (49) 



X=} '^uf^fil), Y 



00 



^?;f*™(J). (50) 



The action I2 remains a constant parameter over a time 
scale of 1/0 (ee'**). 

The angle wi can become a rotating or librating angle 
depending on the initial value of /xj„(/). Stars will be 
near a mean-motion resonance if wi librates. For sim- 
phcity, I introduce the slow angle 9 = Wi -\- — n, 

which casts IC{m, I) into the form 

lC{m,l) ^no{I)-rnnph - ee^*A^(/) cos 6*. (51) 

The distribution function (DF) / = /o + /i is thus con- 
served along the trajectories determined by 



dh 

dlC{m, I) 



dh 



ee"M'„^(/) sinf 



cos 61, (52) 



(53) 



These are equations of perturbed orbits bound to the Ith 
Fourier component in the linear regime. According to the 
dynamical mechanism presented in 21 the / = Fourier 
component is capable of keeping /im(/) ~ O (ee**), which 
in turn, can lead to the resonant capture of stars in the 
linear regime once 9 begins to evolve slowly. 

I have frozen the (small) exponential factor ee''* and 
plotted the instantaneous isocontours (IICs) of /C(2, 0) 
in Figure [4] for the initially circular orbits {I2 = 0) of 
mode S2. The isocontours have been displayed in the 
coordinate plane of ui = sin0 and U2 — \/27i cosf? 
where the transformation {9,Ii) (ui,M2) is canonical. 
Note that the angle 9 is measured clockwise with respect 
to the positive U2-direction. 

The IICs of Figure [4] resemble the topology of the 
celestial three-body problem (e.g.. Wisdom 1980; Win- 
ter & Murray 1997). Two homoclinic loops that inter- 
sect at the instantaneous hyperbolic point ho surround 
a crescent-like, resonant region where the angle 9 li- 
brates. The instantaneous elliptic point co corresponds 
to stable orbits at exact resonance. For a specified I2, 
the coordinates of cq and ho are {9,Ii) = [tt, r7r(t)] and 
[9,Ii) = [O,ro(t)], respectively, where r^(i) and ro{t) are 
the real roots of 



0. stdAl,{r^,l2) 
tJ.,n{rz,l2) - ee 



0, 



0,^, (54) 



at a given time t. 



UNSTABLE DISK GALAXIES. II. 



7 




Fig. 4. — The instantaneous isocontours of the averaged Hamiltonian A2(2,0) for mode S2. The factor ee** is dealt with as a constant 
and has been set to 0.01 {left panel), 0.05 {middle panel) and 0.1 {right panel). The maximum of fCi has been normahzed to unity. The 
crescent- hke resonant zone, driven by the instantaneous eUiptic point cq, has been shaded. Darker regions correspond to smaller values of 
K{2,0). The separatrices intersect at the instantaneous hyperbolic point hg. 



To this end, I show that the eUiptic point Cq always hes 
on the negative U2-axis with 9 ~ tt. I assume the smah 
variations 9 — 9 — n and /i = /i — VT^it) for 1 = 0, and 
hnearize equations ([52|) and ((53|) to obtain 



d9 
'di 
dii 
It 



= ai2{t)Ii, 
--a2i{t)9-r^{t), 



where the time-dependent coefficients are defined as 



ai2{t) ^ 



dh dll 
a2i{t) = ee^' [^^ W] ,,^.^(,) ■ 



r(t) 



(55) 
(56) 

s 

(57) 
(58) 



The quantity rT^{t) is determined by diff'erentiating (|54p 
with respect to t. I obtain 



rn{t)-- 



ai3(0 = 



ai2{t) ' 



dh 



(59) 
(60) 



Eliminating Ii from the linearized equations leads to 

ao{t)9 = a,,{t)+0{e'e^^') 
Mail 



ao{t) 



dh 



(61) 
(62) 



/l=r,(t) 



The sign of the time-dependent spring coefficient ao(t) 
is determined by the sign of dji^/dh- In cored stellar 
disks whose equilibrium potential fields are monotonic 
functions of i?, the following conditions always hold 



dJd, 



(63) 



The equality sign corresponds to the galactic center. The 
isochrone, Kuzmin-Toomre and cored exponential disks 
(investigated in JH and Paper I) fulfill which im- 

plies dfi'^jjdli < 0. Therefore, the spring coefficient 
ao{t) is positive and the homogeneous solution of (|6ip 



is bounded. This proves that the instantaneous critical 
point {9,Ii) = [7r,r7r(i)] is of elliptic type. 

The resonant zone disappears for I ^ when the pat- 
tern speed lies in the interval ([5^ . This is because all 
participating stars in the pattern formation satisfy (|36p 
that yields 



t,'M^li^R + 0(ee^ 



(64) 



Accordingly, /i5„(/) cannot flip sign for I ^ and equa- 
tion (|54p will not have any real root for small pertur- 
bations of ©(ee**). In such a circumstance, IICs will 
constitute a bundle of unidirectional closed curves that 
encircle the origin. In %.'S[ I will verify for the spiral 
mode S2 that pattern stars satisfy ([M)) . 

5.1. The Expansion of Resonant Zones 

It is evident that the structure of IICs in Figure [4] 
evolves with time as the perturbations grow proportional 
to ee**. For small perturbations the critical points cq 
and ho are preserved although they are displaced. As 
a consequence, there always exists a resonant zone for 
ee"* <C 1 and it would be interesting to learn how that 
zone expands in a growing mode. To answer this question 
analytically, I measure the variation of the quantity 

A = do - d„ (65) 

where do and di are, respectively, the values of h on the 
outer and inner homoclinic loops at = tt. The implicit 
functional form of homoclinic loops is given by 

:^o(/i, /2) - ee^*<(/i, /2) COS0 = 

lCo[roit),l2]-ee'*Al[ro{t),l2], 

whose temporal derivative at ^ = tt results in 

dAt 



(66) 



dj. 



dh 
dAf 



+ ese^M^(d„/2) 



dh 



Ii=d„ 

— ese' 

Il=ro 



*<(ro,/2), (67) 



for v = o,i. The bracket on the right hand side of equa- 
tion ([ST]) vanishes because of and the bracket on 
the left hand side is simply d9^/dt {v=o,i). From (|52p 
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and (|53)) one can verify that dOo/dt < and dOi/dt > 0. 
Equation ([67]) thus becomes 



d,, = -- 



[Alid,,l2)+A'i{ro,l2)], y = o,i. (68) 



Subtracting di from do leads to 
'Al,{do,h) 



A = s(ee''*) 



(69) 



The angular rates and 9i can be estimated using the 
linearized equation (|55p as 



^o = ai2(0K-r,(i)], ^» =al2(^)M^-^■^(^)]■ (70) 
Since the resonant zone is thin, one may assume rT,{t) « 
{do + di) /2, which can be combined with ((70|l to obtain 



|ai2(i)|y. 



(71) 



Substituting from ([HI) in ([Ml) yields 



|ai2(t)| 

(d„/2)+^^ {doj2) 



(72) 



where the quotient of the barcketed terms and |ai2(t)| is 
of 0(1). Integrating ((72|) results in 



A 



2Vee'^*C'(l), A « sVee^*C'(l), (73) 

from which the expansion rate of the resonance width 
per unit length is estimated: 



A/A-s/2. (74) 

This shows that stars are steadily captured into reso- 
nance by the I — Fourier component as the mode grows 
in the linear regime. 

My calculations show that the maximum of A^{I) in 
unstable modes does not coincide neither with /ig nor 
with Co, but it falls inside the resonant zone within a 
time scale ^ 1/0(66**) as A grows from zero width. 
The observed shift from the location of cq depends on the 
gradient 9/o/9J0. All growing modes obey the following 
rules near cq: 



oh 



dh 



>o>,lii) if ^>o. 



(75) 
(76) 



In both cases the magnitude of fi'^{I) remains small near 
the maximum of A'^{I), confirming the theory of SJHthat 
the perturbed density profile is supported by the syn- 
chronous precession of orbital axes. 

5.2. Angular Momentum Transfer 

The location of the resonant zone in the (ui, U2)-plane 
provides valuable information of the behavior of resonant 
stars. The Hamiltonian varies according to the equation 

dlC{m,l) dlC{m,l) 



dt 



dt 



f.se''A'^{I) cose. (77) 



One can thus determine how the orbital energy of 
trapped stars changes over time. Recalling that the in- 
stantaneous elliptic point cq occurs at 9 — tt, for most 
resonant orbits the minimum of 9 is larger than 7r/2 (see 
Figured]) and equation ([77]) yields 



dlC{m, 0) 
di 



> 0. 



(78) 



In other words, stars gain energy during a librational 
motion around 9 = 7:. Since I2 = Jr (for ^ = 0) is an 
adiabatic invariant in the resonance zone, an increase in 
the orbital energy of stars boosts their angular momen- 
tmn. 

There is a different mechanism for the angular momen- 
tum transfer to/from stars whose 9 is rotational. Such 
stars will definitely experience both /i > and /i < 
states over a complete period of 9, but which one over- 
whelms the other and how does it affect the orbital an- 
gular momentum? Outside the resonant zone we have 



ImLWI » 



St ^^m(-f) 



dh 



(79) 



which guarantees the rotation of 9 according to ([5^ . I 
now divide equation (|53p by (|52p and ignore all terms of 
0(e2e2''*) to obtain 



dh 
'd9 



mL(^o) 



exp 



s9 



mL(^o) 



Aj,,(/o)sin0, (80) 



where Jq is the initial value of the action vector / = 
(Iijh) a.t 9 = 0. For a clockwise rotation of phase 
space flows, 9 ranges from to 2tt and reversely for their 
counter-clockwise rotation. Moreover, the direction of 
rotation is determined by the sign of iil^{Io). Given 
these points, integrating (|80)) over a complete cycle of 9 
results in the incremental change of h as 



Ahil)^eAi 



sign (/i^ 



/2tts\ 


- 1 









(81) 



It is seen that A/i is positive for > and stars 

gain angular momentum as the perturbations grow. The 
opposite happens for ^5„(Jo) < 0. Combining ((64)) and 
(|8ip shows that the angular momentum gain or loss is 
decided by the sign and magnitude of the Fourier num- 
ber I. The Fourier components with I < and I > 
respectively drain and boost the angular momentum of 
stars so that 

A/i(-l) < A/i(-2) < A/i(-3) < • • • < 0, (82) 
A/i(-f 1) > A/i(+2) > A/i(+3) > • • • > 0. (83) 

Not all stars bound to the / = component are in 
the resonant zone, specially in the limit of ee''* 0. 
Hence, the sign of A/i(0) is determined both by reso- 
nant and non-resonant stars. Based on my arguments 
presented after equation (|77[) , resonant stars always gain 
angular momentum and have a positive contribution to 
A/i(0). However, It is the sign of dfo/dJ^ that deter- 
mines whether non- resonant stars of the I = component 
are emitting or absorbing angular momentum. Equations 
(TTS]) and ([76]) together with (j81|) show that non-resonant 
stars of the I = component gain and lose angular mo- 
mentum for dfo/dJ^ < and dfo/dJ^ > 0, respectively. 
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Fig. 5. — The density components T,\ (top panels) and the contours of ^2(0) (bottom panels) for mode S2. In all figures dashed lines 
mark the stars at exact CR with /i2(f2) = 0. The line /i^^(f2) = lies outside the plot range because it is parallel to the Sl^-axis with 
= 0.454. Straight solid line in the bottom-right panel is defined by ^\ (f2) = 0. Its relative location with respect to the highly populated 
regions of the frequency space shows that stars bound to the I = 1 component evolve far from the OLR. The contour levels of ^2(0) range 
from 10% to 90% of the maximum with increments of 10%. 



Therefore, A/i(0) is certainly positive if the mode de- 
velops in a region of the phase space with dfo/dJcj, < 0. 
Since most stars of the / = component are non-resonant 
as ee''* 0, the occurrence of a growing mode as a result 
of dfo/dJ^ > will lead to A/i(0) < 0. The results of 
the above analysis are consistent with the results of JH 
and the bar charts of Figure [1] 

5.3. Mode Components in the Frequency Space 

To understand how the mechanism of resonant trap- 
ping operates on non-circular orbits in the frequency 
space, I have plotted the components of mode S2 in 
Figure [5l Top panels display the pattern components 
in the configuration space, and in bottom panels I have 
shown the isocontours of ^2(1!) = A2{I) together with 
the lines fi^ii^:) — (straight solid line) and H2{^) ~ 
(dashed lines) . Darker regions in the bottom panels cor- 
respond to larger values of ^2(1!). The contour plots of 
^2(11) show which stars in the frequency/action space 
are engaged with the Ith component. As one could an- 
ticipate, the highly populated regions of all components 
are close to the CR zone of Figure |31 The maxima of 
A2(r2) have not been located on the line H2i^) — be- 
cause the initial density gradient of the equilibrium state 
displaces the center of mass of pattern stars. In modes 
with dfo/dJ^ < near the CR zone, the maxima are 
shifted to regions with /i2(^) > (see ^5.ip . 

Whilst the I — component has trapped near-circular 
orbits into the CR (Figure O, the stars bound to other 
components evolve far from mean-motion resonances 
and the angle 9 becomes rotational for them. In fact, 
stars bound to the I = —1 component avoid the line 
/i,^j^(r2) = and the quantity /x~^(r2) is always negative 
as long as inequality ([37)1 holds. Through a similar mech- 



anism, stars avoid the line /i^(Jl) = and the quantity 
/!,;'„ (H) remains positive regardless of the magnitude of 
ilp. These results show that the dynamical mechanism 
for growing modes of finite s differers with LBK's the- 
ory established in the limit of s — > 0. In other words, in 
a growing spiral mode similar to mode S2 of this paper, 
the I — ~1 and / = 1 components are not associated with 
the ILR and OLR. A new dynamical origin of instabil- 
ities, which interprets the results of this paper and has 
not the restrictions of LBK's theory, is presented in the 
next section. 

6. THE ORIGIN OF INSTABILITIES 

Consider a randomly generated, small-amplitude, m- 
fold density wave of pattern speed Qp whose potential 
field can be generally expanded in the Fourier series 
of the angle variables. Independent of the form of ini- 
tial phase space distribution, determined by /o(«/), the 
points (1,9) = (/o,0) and (1,9) — {Io,t^) will emerge 
as time- invariant stationary points of the flows generated 
by /C(m,0) if the following condition holds 



dh 



(84) 



at / = In such a circumstance, orbits associated 
with either of these stationary points will never change 
their action vector / — (/i,/2) and their orbital energy 
must be conserved. This condition will be satisfi e d only 
if s = in 173. This is how a Ivan Kampen I (jlQSSl ) 



mode is born. The state of (|84p can occur at any point 
in the infinite dimensional action space and for arbitrary 
values of ilp. Consequently, stationary modes constitute 
a continuous family and Mathur's (1990) isolated, pure 
oscillatory modes are not feasible in stellar disks. 
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Let me now suppose that for some pattern speed fip 
in the interval ([55]) . the quantity n'^i^) accidentally re- 
mains small (for a group of stars) over a finite duration 
of time, but the condition (|84p is violated. This (likely) 
symmetry-breaking phenomenon creates a slim resonant 
zone of the width A(i) (see There is not any phys- 

ical constraint on the evolution of A(t) when the pertur- 
bations are in their early stages. For a steady rate given 
by (|74p more stars are trapped by the resonant zone if 
s > and the density wave corresponding to the Fourier 
number / = is magnified. As I discussed in ^5.2l the 
angular momentum of the captured stars increases and 
some other reacting stars should therefore respond in or- 
der to recover the angular momentum balance of the disk. 
The reacting stars live in the same CR zone defined by 

As more stars are trapped into resonance by the I = 
component, non-resonant components with |Z| > should 
also involve more stars to compensate the angular mo- 
mentum deficit. So the amplitudes of other wave com- 
ponents increase as well. This is what I am considering 
as the origin of instabilities: triggering unsteady den- 
sity perturbations by the capture of stars into the CR. 
Angular momentum transfer between non-resonant stars 
cannot be facilitated if a resonant gap does not open in 
the phase space. In fact, the initial destabilizing imbal- 
ance of £ is generated by an irreversible engagement of 
resonant stars even if the amount of angular momentum 
that they absorb is small (e.g.. Figure 8 in JH). A per- 
turbation with the property of ([5^ fixes s — 0, does not 
awake stars of |Z| > components, and it creates a sta- 
tionary mode. I conclude that most (if not all) unstable 
modes disappear for fip > fig just because the CR is 
destroyed in such a circumstance. 

The I — component is not always an angular mo- 
mentum absorber. For example, in models with an inner 
cutout of the DF, one has dfo/dJcf, > which enforces 
Lm{0) < (discussions in §5.21 and JH). Figures 2 and 
11 in JH demonstrate two modes that show angular mo- 
mentum emission by the I — component. The pattern 
speeds of some C-modes explored in Paper I exceed flo 
marginally. Those modes also have Lm{0) < 0. Unstable 
modes with flp > VIq are barely observed in modal cal- 
culations because their large pattern speeds shrink them 
to the galactic center. Coexistence of mode Bl and a 
compact C-mode can explain the origin of double-barred 
galaxies. For a softened-gravity model of the exponential 
disk with an inner cutout. Alar Toomre (private commu- 
nication) also finds an inner edge mode that has a corota- 
tion circle inside its pattern. Toomre's inner edge mode 
is another example of a mode with df^/dJ^ > and 
L m jO) < 0. Axisymmetric features /grooves introduced 
bv ISellwood fc Kahn I (|1991[ ) to the surface density and 
DF also develop a spiral mode whose pattern is discon- 
tinuous at the location of the groove. That discontinuity 
is due to the rapid change in the sign of dfo/dJ^: while 
the I = component boosts the angular momentum of 
stars outside the groove, it has an opposite effect on inner 
stars. 

Each unstable mode occupies a finite region in the fre- 
quency space due to the finiteness of its CR zone. This 
is how divisions (gaps) are created between the pattern 
speeds of unstable modes, the Jl-space is quantized by 
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Fig. 6. — The resonant zones of modes Bl, SI, S3 and S6 for the 
cored exponential disk with {N,\,a) = (6, 1,0.42). Contour plots 
show the IICs of K{2, 0) for I2 = and ee"' = 0.01. Darker regions 
correspond to larger values of /C(2, 0). Contour levels of each zone 
have been normalized to the maximum value of /C(2, 0) in the same 
zone. 

resonances, and finally, unstable modes constitute a dis- 
crete family. Growth rates are sorted according to the 
likelihood of resonant capture, which is proportional to 
the initial population of stars determined by the equilib- 
rium DF. In a stellar disk with a falling density profile, 
the growth rates of linear modes are sorted in a decreas- 
ing order from compact modes to more extensive ones. 
Most barred and spiral modes obey this rule (see Pa- 
per I) . The perturbation theory of [JS] says that stars live 
on the integral manifolds of constant I2 before and af- 
ter capture into resonance. In order to guarantee a sus- 
tained growth of the CR zone, the initial DF fo{J) must 
be smooth and non-zero along the curves of constant I2 , 
and has to flip its sign there. A resonant zone 

cannot grow in a region of the phase space which is de- 
void of stars from the beginning, or its stars have been 
forced to migrate by unsteady processes. The first-order 
perturbation theory is valid as long as the CR zone of 
each mode is inaccessible by the stars of other modes 
and I2 is an adiabatic invariant. 

7. MODE SATURATION AND STABLE DISKS 

Perhaps the most important issue related to the growth 
of instabilities is that the amplitudes of unstable modes 
are saturated after several pattern rotations. Saturation 
is indeed a nonlinear phenomenon an d can have com- 
plex routes. ISellwood fc Binnevi (|2002D suggest that the 
growth of a mode is stopped by the emergence of horse- 
shoe orbits. What they call "horseshoe orbits" are sim- 
ply the librational orbits captured by the CR. I found 
those orbits the main caus e of instabilities. I agree with 
ISellwood &: Binnevi ()2002l ) that the linear perturbation 
theory must fail at some stage, but it does not necessarily 
imply that resonant zones can self-control their expan- 
sion beyond the linear regime. Using the phase space ge- 
ometry and the nonexistence of the first integral I2 at the 
boundaries of neighboring resonances, I argue that reso- 
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pansion direction of the / = Fourier component as the 
modes grow. The growth rates of modes are different and 
mode S6, which is at the bifurcation point of the S-family 
(see Figure 2a in Paper I), has the smaUest s. The con- 
figuration of the CR zones in the phase/frequency space 
and their thickening according to equation (|74p suggest 
that resonance overlapping is likely in stellar disks with 
rich branches of unstable modes in their eigenspectra. 

To find out how and when an overlap takes place, I as- 
sume two growing waves with the eigenvalues mflp^i+isi 
and mr2p_2 + *S2 that successively appear in the spec- 
trum with r^p^i > i7p.2. These waves develop CR zones 
of widths Ai(/2) and A2(/2) on each manifold of con- 
stant I2, and I define (^12(^2) as the distance between the 
instantaneous elliptic points of those zones. Let me sup- 
pose that the global minimum of di2{l2) corresponds to 
the critical action I2 = ■ An overlap thus occurs if 

dl2(/2")^i[Al(/2")+A2(/2"-)]- (85) 



Fig. 7. — Left panel: Contour plots of A^(f2) for modes Bl, SI, S3 and S6 of the cored exponential disk with {N, A, a) = (6, 1, 0.42). The 
levels of isocontours have been normalized to the maximum amplitude. Right panel: The emergence of ILR through the intersection of the 
lines /Jm^(f2) = with the frequency space. The line /i5^(f2) = (it is independent of m) has also been plotted to show the location of the 
CR. In all cases the pattern speed has been set to Hp = 0.1. There is no ILR associated with m = 1 perturbations for any Qp > 0. 

nance overlapping is an alternative and efhcient mecha- 
nism to stop the growth of unstable modes. The stabi- 
lization of stellar disks in the presence of the ILRs is also 
explained by the same mechanism. I then present a sim- 
ple criterion for the stability of non-axisymmetric pertur- 
bations, and discuss about the stability of self-consistent 
scale-free disks. 

7.1. Resonance Overlapping 

As a mode grows, the width of its CR zone increases 
proportional to V ee"* until two adjacent zones, corre- 
sponding to modes of different pattern speeds, overlap 
and compete for trapping the stars in the overlapping 
zone. This is the moment that a chaotic layer occurs in 
the phase space according to Chirikov's (1979) theory, 
and stars living in that layer repeatedly migrate between 
different resonant zones. Thus, the share contributed to 
each competing resonant zone from stars in the overlap- 
ping region drops substantially and the growth is sup- 
pressed. A natural consequence of the emergence of a 
chaotic layer is the radial migration of stars that heats up 
the disk. Some modes may be gradually dissolved as their 
Kolmogorov-Arnold-Moser (KAM) tori around stable pe- 
riodic orbits are destroyed by modes with larger resonant 
zones and growth rates. I remark that the weighted resid- 
ual form of the CBE that I obtained in Paper I is valid as 
long as we are allowed to average out angle variables from 
dynamical equations. After the occurrence of chaotic or- 
bits it is impossible to average out resonant angles and 
extra coupling terms between the amplitudes of unstable 
modes are contributed to the reduced CBE. Details of 
my numerical simulations in the nonlinear regime and in 
the presence of chaotic orbits will be presented in Paper 
III. 

Figure [6] shows the I — resonant zones of modes Bl, 
SI, S3, and S6 for the initially circular orbits {I2 — 0) of 
the model introduced in S}2l This is the first demostar- 
tion of its kind that shows how the resonant zones of 
different unstable modes partition the phase space of a 
stellar disk. Since the perturbations grow according to 
an exponential law, resonant zones had been slim rings 
infinitely long ago. The structure shown in Figure [6] has 
been obtained by setting ee** = 0.01 for all modes. Cor- 
respondingly, resonant gaps are of comparable size at the 
displayed moment but they will not remain so because 
the growth rates are different. Left panel in Figure [7] 
shows the amplitude function ^^(O) of the same modes 
of Figure [5] but this time in the frequency space that 
covers all possible values of l2- Arrows indicate the ex- 



On the other hand, the magnitude of di2{l2^) can be 
estimated as 



(86) 



di2{If)~ha-h^+0{ee''), 

where = flpj {j—1,2). By combining (185 

and and using ([75]) . I obtain 



< 



Ve^Oi ( 1 ) -f Ve^02 ( 1 ) 



l2=Io 



(87) 



which will be satisfied in the linear regime at some t = t^r 
for a sufficiently small value of /1.2 — Ii.i- If the post- 
linear changes in the DF allow for a sustained expan- 
sion of resonant zones (when the amplitude growth is no 
longer exponential), one can still anticipate an overlap 
for larger values of Ii^2 — Cooling a stellar disk in- 
creases the number of unstable modes (Figure 6 in Paper 
I) and decreases the gaps between the CR zones. Conse- 
quently, transition to chaos is faster in cold disks than hot 
ones. This result is in harmony with Sellwood's (2007, 
private communication) iV-body experiments. The very 
first overlapping does not necessarily occur between the 
modes of the same azimuthal wavenumber to. The pat- 
tern speeds of two rapidly growing modes of different to 
may result in a very small Ii^2 ~ -^1.1 and thus cause an 
early-stage overlap. 

The fundamental bar mode may still be saturated in 
the absence of spiral modes and in a shorter time scale 
than chaotic diffusion of stars. When the resonant zone 
of mode Bl grows, the inner homoclinic loop rapidly 
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shrinks to central regions, the supply of stars to the CR 
zone is cut from the galactic center, and the bar mode 
saturates there. On the other hand, the bar continues 
to elongate because of the expansion of the outer homo- 
clinic loop. The elongation rate, however, will slow down 
if dfo/dJif, < 0. A declining rate of amplitude growth 
gives time to the trapped stars to fill their invariant tori 
and maintain the bar in outer regions. 

7.2. The Role of Inner Lindblad Resonances 

Disappearance of unstable modes when the pattern 
speed violates ([57| can also be explained by the same 
mechanism that saturates instabilities. The ILR emerges 
when the line ii^{n) = crosses the circular orbit 
boundary = rc(f2i?). This happens while the line 
/x^(r2) = has a crossing too and the CR zone has al- 
ready been created. The competition between the ILR 
and CR near their neighboring boundaries can therefore 
prohibit both resonances from broadening their zones 
and the disk develops a stationary wave. Right panel 
in Figure [7] shows the frequency space of the cored ex- 
ponential disk together with the lines — and 
/i^j^(r2) = that have been drawn for f7p — 0.1 and 

m =1,2,3. The line /ij~^(r2) = never intersects the fre- 
quency space independent of the magnitude of r2p > 0. 
A crossing is possible for m > 2 and the slope of 
/^^^(Jl) = increases proportional to m. Subsequently, 
the distance between the CR and ILR zones decreases by 
increasing m and the disk shows more reluctance against 
developing an unstable mode (see Figure 1 in Paper I). 
iBertin et al. I jl977 ) had also reached a similar result for 
their three- and four-armed modes, and interpreted it 
based on the closeness of ILR and CR but with no men- 
tion of the increased likelihood of resonance overlapping. 

Having resonant zones of finite width for both the CR 
and ILR is an essential prerequisite for their subsequent 
competition. The lih Fourier component in the radial an- 
gle can trap stars into a growing resonant zone if /o( J) is 
smooth and non-zero on a range of manifolds of constant 
I2 — Jr. — IJtp/m, and /ij,j(r2) can switch its sign in that 
neighborhood. Having some population of stars on circu- 
lar orbits always guarantees the creation and growth of 
the CR zone because I2 is identical to Jr and Ii = J,j,/m 
can vary independently. However, if the ILR exists, its 
adiabatic invariant is defined by I2 — Jr + J4, /m and the 
variation of Ii in the vicinity of a resonant value implies 
the variation of Jr. Therefore, the ILR can develop a 
growing resonant zone if fo{J) has a functional depen- 
dence on both actions. Stellar disks with such a prop- 
erty will have non-zero radial velocity dispersions (warm 
disks). 

Let me now consider a group of stars moving on near- 
circular orbits that satisfy //"^(fi) > and have been 
suddenly trapped into resonance by the ILR. On the res- 
onant tori, the orbits of those stars are deformed so that 
I2 ^ Jr -'r J(t,/m remains constant. This means that or- 
bits are elongated as Jr continuously increases from its 
tiny initial value and drops until /x^^(r2) becomes neg- 
ative and phase space flows reach their turning points. 
At this moment, some stars are likely to find themselves 
on the integral manifolds of constant Jr^ Q and join the 
CR zone while they carry energy and angular momentum 
away from the ILR. A reverse migration is also possible. 
The switchings can only take place in the overlapping 



region of the ILR and CR where highly eccentric orbits 
live and ^"^(Jl) < < /^"^(Jl). An effective early-stage 
overlap requires a non-zero fo{J) hi the cross section 
n iS~^ of the sets 

s^„,^{n■.^,l{n)^o{e)}, (88) 
s-,'^{n■.^,-^\n)^o{e)}. (89) 

As a result, infinitesimal wave amplitudes (in the linear 
regime) can suppress their own growth due to the compe- 
tition between the CR and ILR. In the overlapping region 
5^n5~^ neither I2 = Jr (for I = 0) nor I2 = jR + J^/m 
(for / = — 1) are conserved and that region is filled by 
chaotic orbits. The chaotic layer should be very narrow 
in the linear regime, but it will bec ome thick in steady 
large-ampfitude bars (jContopoulos 1 11983). 

The stability of warm stellar disks for Qp < r2iLR(m) 
can be deduced from an independent mathematical rea- 
soning as follows. Equation (j74p applies to any growing 
resonant gap in the phase space. Let me assume that 
the ILR and CR of a wave with the pattern speed flp are 
widening their resonant zones with the rates 

Acr/Acr - si/2, Ailr/Ailr - S2/2. (90) 

In an equilibrium disk with cEo/dR < 0, the probabil- 
ity of capture into the ILR is more than the CR, and 
according to (|90)) . we need S2 > si to assure a faster ex- 
pansion rate of the ILR zone. This inequality contradicts 
the basic assumption of linear analysis that separates the 
time-dependent part of perturbations as exp{—iu)t) and 
assigns the same growing envelope and pattern speed to 
all parts of a normal mode. Therefore, the only possibil- 
ity for a disk with a monotonic density gradient will be 
si = S2 = 0. In cold disks with /o(J) = S{JR)go{J4>), 
where S{Jr) is Dirac's delta function, the ILR cannot 
open a resonant gap because it suffers from the lack of 
fuel (stars) needed for a possible expansion along the 
curves of constant I2 = Jr + Jip/m. This automatically 
implies S2 = and a unique value is assigned to si. The 
CR zone can thus grow freely and destabilize the disk. 

It is worth noting that with an inner cutout, the sur- 
face density rises in central regions, reaches to a maxi- 
mum and then decays outwards. Denote i?cR and -Rilr 
as the corotation and inner Lindblad radii, respectively. 
If the cutout is not sharp and the density peak lies fairly 
between -Rilr and i?cR (see the leftmost panel in Fig- 
ure 8 of Evans & Read 1998b), it will be possible to 
find two points of the same density on the curve of 
So(-R): one near the ILR on the rising part of T,q{R) 
and the other near the CR on the falling side, i.e., 
So (-Rilr) ~ So(-Rcr). Thus, the ILR and CR can man- 
age to grow with the same rate si = S2 and an unstable 
mode emer ges. Such modes have b een demonstarted in 
Figure 7 of lEvans fc Read I ()1998bD whose cutout func- 
tions immobilize the circular orbits of the galactic center 
and all near-radial orbits. Accordingly, the population 
of active stars in the overlapping region of the ILR and 
CR declines substantially and the simultaneous growth of 
both resonances becomes likely. By decreasing the pat- 
tern speed, the ILR moves outside the cutout radius and 
its resonant zone becomes visible to the CR. The stellar 
disk is therefore stabilized. The condition si = S2 7^ 
cannot be reached if So(-Rcr) differs remarkably from 
So (-Rilr)- This happens when both the ILR and CR lie 
outside the cutout radius or when the cutout is so sharp. 
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The above analysis shows how diverse the response of 
disk galaxies m ay be to non-axisymmetric excitations. 
iToomre I ()1981|) suggested that the two parameters Q and 
X decide the fate of non-axisymmetric density waves. 
Wave mechanics based on the behavior of resonances, 
shows more complex decisive factors that govern the evo- 
lution of density waves, specially near the critical values 
of Qp « J7ilr(''ti) and for non-monotonic gradients of 
equilibrium density. 

7.3. A Criterion for Global Stability 

I am now in a position to provide a simple criterion for 
the stability of stellar disks in the absence of the ILR. 
It is the CR that captures stars, contributes an angular 
momentum imbalance to the stellar disk, and triggers 
unstable modes. Circumferential waves will be globally 
stable if resonant gaps do not open near the CR for all 
possible pattern speeds. Consequently, Lm{0) should re- 
main zero, which according to equation (j2ip implies 



0, yj,k>o. 



(91) 



This condition is satisfied only if [see equations ((T6)l and 

m] 

^ = 0, for aU J e 3?+ X 3?+. (92) 

OJd, 



The stability criterion given in ([92)1 can also be deduced 
from the earlier works of LBK and JH as I explain below. 
An alternative form of ([^0]) is given in equation (B6) of 
JH as 



dt ^ 

1 — — OC 



dj I 



dfo , dfo 



dJd, 



where 



V5 ^^a^^f (J). 



(93) 
(94) 



Dividing equation (|93p by 167r^ and changing the sign of 
u) lead to equation (28) in LBK. The amount of angular 
momentum emitted/absrobed by the I = component is 
therefore determined by 



d£^ 



= -2{rrmfse^ 



1=0 



dJ 



dfo 



\Vn\ 



dJrj) \mfl^ 



■to 



(95) 



which reduces to equation (96) of lGoldreich fc Tremaine 1 
()1979| ) as s ^ 0. It is seen that the sign of dfo/dJ^ de- 
termines whether the angular momentum is added to, or 
drained from the I — wave component. Furthermore, 
the CR will not imbalance £ if the condition ([92| holds. 
This means that the condition (IM)) cannot be violated 
according to my discussion of fJH and s should vanish for 
arbitrary perturbations. A stellar disk that satisfies (|92p 
will be hot because it is impossible to reproduce the sur- 
face density gradient in the limit of — > with a DF 
that does not depend on J^. It has now become more 
transparent why there is a correlation between Toomre's 
(1964) Q and the disk stability for non-axisymmetric ex- 
citations. An illustrative example has been given in Fig- 
ure 6 of Paper I. Most of S-modes disappear when the 
parameter N of the DF is decreased and more stars are 
distributed on near-radial orbits. 



7.4. The Stability of Scale-Free Disks 

Stability analysis of cuspy ma ss mod e ls has been a big 
puzzle in stellar dynamics since [Zang I (|1976f ) computed 
the global modes of Mestel's (1963) disk and showed the 
importance of to = 1 excitations. To tackle the sin- 
gular center of Mestel's disk, he introduced an inner 
cutout that freezes central s tars w ith diverging orbital 
frequencies. 'Evans & Read (1998a,b) followed the same 
procedure but for different cusp slopes. In spite of all 
these remarkable efforts, what H02 illustrated for the fre- 
quency space of scale-free potentials was indeed the key 
to resolve the long-standing stability problem of scale- 
free disks without any simplifying cutouts. The map of 
the angular sector that H02 plotted in his Figure 1 for 
the scale-free potentials (jTouma fc Tremaine Ill997[ ) 



VoiR) 



sign(/3)i?'3, /3^0, 
Ini?, /3 = 0, 



1< /3 < 2, 



(96) 



will be an angular sector in the (r^ij, r2i)-space with the 
ri/j-axis and the straight line 



1 



1 



(97) 



being its lower and upper boundaries, respectively. The 
slope of the upper boundary is always less than or equal 
to 1/2. It is therefore obvious that except for to = 1, all 
/i,^j^(r2) = and /i°„(r2) — lines will intersect the fre- 
quency space (whatever ftp may be) and resonance over- 
lapping between the CR and ILR can occur if there are 
enough number of non-circular orbits. In other words, 
the radial velocity dispersion du (defined in Evans & 
Read 1998a) should exceed a critical value (7u,min{m) to 
give rise to nonempty cross sections n S~^. This 
is possible only for to > 1 excitations and it suggests 
that scale-free disks are subject to to = 1 instabilities. 
Finding au,mini'm) is an open, yet conceivable, problem 
that can generalize Toomre's (1964) criterion to non- 
axisymmetric perturbations of scale-free disks. I ex- 
plained in §7.2l that the ILR cannot help to stabilize cold 
disks. That argument applies to self-consistent scale- 
free disks as well, and strongly justifies the existence of 
o'M,min('7i) for TO > 1. Morcover, a stability condition 
like du > o'„,min(?Ti) may correlate with Lynden-Bell's 
(1993, equation 6-23) criterion, which has been formu- 
lated based on the angular dispersion of lobe tumble 
rates. 

8. DISCUSSION AND CONCLUSIONS 



I had already given an evidence in gSTSlthat LBK's the- 
ory cannot be extended to large growth rates. Here I pro- 
vide a more detailed analysis of LBK's mode mechanism. 
Taking the s ^ limit of (|93p leads to equation (30) of 
LBK. The resulting integrand of the Zth component in- 
volves the term S{lflji + mfl(i, — u!) but its argument never 
becomes zero for ^ < — 1 as long as flp exceeds f2iLR(TO). 
Once the delta function vanishes for Qp > r2iLR(TO), the 
angular momentum content of ^ < —1 components re- 
mains zero. Therefore, angular momentum transfer be- 
tween inner and outer resonances as suggested by LBK, 
is feasible only for Hp < Qii,^ {m). This is a constraint 
on the applicability of LBK's theory. 

According to the arguments of H7.'2\ and the WKB 
theory, the ILR is not transparent to short wavelength 
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disturbances with X ^ 1 (|Binnev fc Tremaine 1 120081 : 
iMark I [l97l and any developing spiral structure must 
be damped unless the ILR remains invisible to the CR. 
Therefore, an inconsistent point in LBK's paper is the 
imagination of an angular-momcntum-transferring spiral 
structure while the condition Qp < f2iLR(?n) has already 
abandoned the existence of such structures. An issue yet 
needs to be explained: If there is no spiral structure in 
stable disks, which mechanism does transfer the angu- 
lar momentum between inner and outer resonances? In 
fact, in the limit of s — ^ the density components 
and S']"*" become stationary waves that are azimuthally 
separated by a phase shift of 90° (like the components of 
mode Bl in Figure [2]). Thus, they exert opposite grav- 
itational torques on each other without any need for a 
communicating spiral structure. 

In this paper I decomposed unstable modes of a model 
disk galaxy to its constituent Fourier components and 
showed how different components experience a gravita- 
tional torque. My results clearly showed that only the 
Fourier component associated with the CR generates a 
resonant zone in the phase space, and other components 
exchange angular momentum far from resonances. This 
result led me to introduce a new dynamical mechanism 
that triggers unstable modes. According to my calcula- 
tions, an irreversible resonant capture of stars into the 
CR causes a synchronous precession of their orbital axes, 
which in turn, support a rotating density pattern. The 
emerged pattern grows because the resonant zone ex- 
pands in the frequency space. 

The irreversibility of resonant trapping is the most de- 
structive event that happens in a stellar disk when a 
group of resonant stars with di < h < do gain angular 
momentum forever. An immediate consequence of this 
phenomenon is the angular momentum imbalance in the 
disk, which requires proper reaction of other stars in or- 
der to conserve C. This is how the Fourier components of 
/ 7^ are generated and the angular momentum transfer 
between them is initiated. Reacting stars do not have a 
librating 9 and they cooperate in a way that the quantity 
remains small for all components. 

Using the maps of resonance zones in the phase and 
frequency spaces, I argued that resonance overlapping 
should stop the growth of unstable modes and stabilize 
stellar disks in the presence of an ILR. I showed that for 
m > 1 the emergence of an ILR is inevitable in scale- 
free models. The competition between the CR and ILR 
can therefore stabilize suffi ciently warm sc ale-free disks 
against m > 1 excitations. lEvans fc Read I ([l998b) sug- 
gested in §6.1 of their paper that self-consistent scale- 
free disks (without inner cutouts) do not admit growing 
non-axisymmetric modes at all, and they ruled out the 
possibility of a critical temperature. Their prediction is 
in agreement with my results except in two aspects: (i) 
Since the line = does not intersect the fre- 
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quency space of scale-free models, there is no ILR to 
compete with a growing CR and density waves will be 
amplified for m — 1. (ii) The ILR in cold disks is very 
special and it cannot develop a resonant zone of finite 
size. It is therefore hard to believe a serious influence by 
the ILR on the phase space flows up to a critical temper- 
ature. 

Note that each unstable mode dominates an isolated 
region of the frequency (action) space until a chaotic 
layer emerges due to resonance overlapping. The rate by 
which the chaotic layer diffuses itself in the phase space, 
is determined by the resistance of KAM tori against com- 
peting resonant zones. The limited space that the fastest 
growing bar mode occupies after its saturation (Khop- 
erskov et al. 2007), is an A^-body evidence for such a 
resistance of spiral modes that fill outer regions of the 
cored exponential disk. 

The other implication of resonance overlapping, which 
has observational support too, is that the central bar in 
grand-design barred spiral galaxies must join the spiral 
arms at the tips of the bar, which are the overlapping re- 
gions of two neighboring B and S modes in the frequency 
space (see Figure [7|). I note that the pattern speeds of 
the spiral and bar components are not the same because 
these structures are associated with different modes. The 
overlapping region is filled by chaotic orbits that yield 
the turbulent mixing of density waves. This process can 
enhance star formation near the tips of the bar. If we 
accept this scenario, barred spiral galaxies should have 
been formed from hot stellar disks whose eigenfrequency 
spectra include few spiral modes. On the other hand, 
flocculent spirals with many wave packets will be the 
natural destiny of initially cold disks that give birth to a 
rich family of spiral modes. 

The mathematical background for the nonlinear evolu- 
tion of modes was developed in Paper I but that formu- 
lation is valid as long as averaging is allowed over angle 
variables. In the presence of chaotic orbits it is impos- 
sible to average out resonant angles and some modifica- 
tions are required to deal with the CBE in its full nonlin- 
ear form. I will present such modifications in Paper III 
with the aim of discovering the dynamical processes that 
generate different classes of barred and spiral structures. 
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